Not to appear in Nonlearned J., 45. 



Relativistic particle acceleration in developing Alfven turbulence 

S. Matsukiyo and T. Hada 

Department of Earth System Science and Technology, Kyushu University, 
6-1 Kasuga-Koen, Kasuga, 816-8580, Fukuoka, Japan 

matsukiySesst . kyushu-u .ac.jp 
ABSTRACT 

A new particle acceleration process in a developing Alfven turbulence in the 
course of successive parametric instabilities of a relativistic pair plasma is in- 
vestigated by utilyzing one-dimensional electromagnetic full particle code. Co- 
herent wave-particle interactions result in efficient particle acceleration leading 
to a power-law like energy distribution function. In the simulation high energy 
particles having large relativistic masses are preferentially accelerated as the tur- 
bulence spectrum evolves in time. Main acceleration mechanism is simultaneous 
relativistic resonance between a particle and two different waves. An analytical 
expression of maximum attainable energy in such wave-particle interactions is 
derived. 

Subject headings: particle acceleration, parametric instability, developing turbu- 
lence, relativistic effect 



1. Introduction 

Large amplitude Alfven waves are ubiquitous in space and astrophysical environments. 
Such waves are beheved to play key roles in the so-called DSA (diffusive shock accelera- 
tion) process in which charged particles are diffusively accelerated in the course of multiple 
scattering through turbulent Alfven waves upstream and downstream of coUisionless shocks. 
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implicitly assumed that the Alfven turbulence is phase randomized and its spectrum is time 
stationary. On the other hand, turbulent Alfven waves commonly observed in the solar- 
terrestrial environments are often intermittent, and coherent MHD structures are frequently 



superposed. Naniely, the phase random a pproximation cannot be assumed (Ide Wit et al. 



I999I : iHada et al.ll2003l: 



Narita et al 



20061) and a spectrum of turbulence may evolve both 



in space and time (IBruno fc Carbondl2005l ). This is probably due to the fact that nonlin- 
ear wave-wave interactions tend to generate coherence among wave phases, and that spatial 
and temporal scales of relaxation processes in space plasmas are much larger than typical 
scales of our solar-terrestrial system. This may hold true with some high energy astro- 
physical environments. That is, spatial and temporal scales of a relaxation process are not 
negligibly small in comparison with scales of a whole acceleration site. Generally speaking, 
particle acceleration rate through coherent wave-particle interactions i s much higher than 
that throu gh incoherent, or fully turbulent, wave-particle interactions (iKuramitsu fc Hada 
2OOOI . I2OO8I ). Therefore, acceleration processes through coherent wave-particle interactions in 
a 'developing' turbulence, where a turbulence has not been fully developed, should be paid 
more serious attention. 

It is well-known that in various space plasma environments a nonequilibrium ion dis- 
tribution function generates large amplitude Alfven waves through some instabilities and 
that those waves nonlinearly evolve to produce coherent magnetic wave forms. One of 
the common interpretations of generation mechanisms of such wave forms accompanying 



density fluctuations whic h are frequently observed in the solar wind (ISpangler et al.lll997 



Spangler &: Fuselierlll988l ) is a parametric instability where nonlinear wave-wave interactions 



convert energ y of a parent wave into several daughter waves wit h different frequencies and 



wavelengths (iGaleev fc Oraevskiil Il963l : ISagdeev fc GaleevI 119691 ) . Although in the context 



of cosmic ray acceleration in turbulent media the above process was taken into account 
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teractions have never been considered. Recent numerical studies revealed that successive 
parametric instabilities re s ult in turbulent wave forms w hich have not been fully developed 
jMatsukivo fc Hadal I2OO3I : iNarivuki fc Hadal boosl . bood l However, only a few past stud- 



ies paid much attention to particle acceleration processes in such a developing turbulence. 
Furthermore, only limited spatial and temporal evolutions of wave forms or spectra were 
discussed, since a computational resource was limited. 

In this paper long time evolution of parametric instabilities of a large amplitude Alfven 
wave in a rather large spatial domain is reproduced by utilizing one- dimensional relativis- 
tic full particle-in-cell (PIC) code. A plasma is assumed to be composed of electrons and 
positrons, since electron-positron pairs can be the dominant constituent of some high en- 
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ergy astrophysical plasmas like in the vicinity of a pulsar, active galactic nuclei (AGN), 
gamma ray burst (GRB), and so on. In the simulation we observe a particle accelera- 
tion process which is quite efficient and is due to interactions between coherent Alfven 
waves and relativistic particles. The process is quite different from some other coherent ac- 
celeration processes being discussed recently, in which high frequency electrostatic waves 
play essenti al roles, like electron sur fi ng acceleration induced by a cross field Buneman 



instabilitv (IShimada fc Hos 



Dieckmann et al.ll200 
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19831 : lLvubarskvl[2006l : iHoshinol [20081 ) 



2002 



Amano fc Hoshinol 20071) and wakefiel d acceleration ( Tajima fc Dawson 



Simulation settings and results are represented in section 2. The acceleration process 
is discussed in detail by using an analytical model and test particle simulation in section 3. 
Summary and discussions are given in section 4. 



2. ID PIC Simulation 

Long time evolution of parametric instabilities of a large amplitude monochromatic 
Alfven wave in a relativistic pair plasma is reproduced by performing a one-dimensional PIC 
simulation. A parent wave is given only at the beginning of the run as a monochromatic and 
right-handed circularly polarized Alfven wave with amplitude BjJBq = 1 and wavenumber 
koc/Qo = 6.21 (number of wave mode is 512) which gives a frequency ujq/^Iq = 0.627, 
where Bq is the ambient magnetic field which is along the x— axis, c denotes speed of light, 
and Qq = eBo/mQC is nonrelativistic gyrofrequency, respectively. Correspond ing velocity 



pertu rbations of the pair plasma are given by the relativistic Walen relation (IHada et al. 



2004j ). Boundary conditions are periodic for particles and all field components. System size is 
L = 518.1c/f2o. Squared ratio of nonrelativistic gyro and plasma frequencies is ft^/ujp = 0.1, 
and normalized scalar temperature for both electrons and positrons is T /ttiqc^ = 1.6 x 10^"^. 

FigH] shows time evolutions of energy densities (upper panel) and Fourier spectra of 
wave amplitudes for and field components (lower panels). In the upper panel the 
solid lines denote electron and positron kinetic energies (labeled by 'elec' and 'posi.'), 
transverse magnetic field energy ('i?^'), and transverse and longitudinal electric field energies 
{^E±^ and '-Ex'); respectively. The dashed line indicates the total energy, which is well 
conserved during the run. Although rather long system evolution (up to QqI = 1616) is 
calculated, the system is still far from the so-called equilibrium state. In the lower panels 
wave amplitude spectra for positive and negative helicity (corresponding to positive and 
negative wavenumbers respectively) modes of the B^ component are shown in the left (B?) 



and middle (-B^ ) panels by using a technique of Fourier decomposition (iTerasawa et al. 
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19861 ). Note that in the middle panel the actual sign of the wavenumber is negative. Since 



most of the daughter waves are right-hand polarized, the left (middle) panel shows wave 
intensity mainly of positively (negatively) propagating waves. At the beginning only the 
parent wave has significant intensity at kc/Qo = 6.21 in the left panel, although it is not 
particularly outstanding in the figure because of the narrow spectrum. In the early stage, 
50 < flot < 100, the parent wave energy is transfered to a variety of daughter waves through 
large number of channels of wave- wave interactions. In this stage two parametric instabilities 
are dominant. A modulational instability generates daughter magnetic fluctuations with 
the same signs of wavenumber as that of the parent wave. They are observed in the 
spectrum. Couplings among the daughter and the parent waves result in compressional 
electrostatic fluctuations with wavenumber smaller than that of the parent wave as seen in 
the right panel. In a decay instability, on the other hand, daughter magnetic fluctuations 
have negative wavenumbers which appear in the spectrum, while daughter electrostatic 
fluctuations have wavenumbers larger than that of the parent wave in the right panel. Most of 
the parent wave energy is wasted in this stage. However, the intensities of the daughter waves 
are still large so that successive parametric instabilities occur thereafter. These successive 
processes are sustained mainly by decay instabilities which can be confirmed because peak 
wavenumbers of the Ex spectr um is always larger than th ose of the B^ and B^ spectra. The 
previous simulation study by iMatsukiyo fc Hadal (120031 ) confirmed occurence of up to the 
second decay instability for the case of Bp/ Bq = 0.1. 

FigJ2^ demonstrates electron energy distribution functions at various times. Fourier 
amplitude spectra of B^ component at the corresponding times are plotted in FigJ2)D. Rapid 
heating occurs after the first instabilities developed {fiot = 95: black dashed line). Up to this 
stage, most of the electrons show semi-stochastic motions in a noizy system with a primary 
wave. Here we refer the primary wave and the noize to the parent wave and superposition 
of the daughter waves, respectively. This results in rapid scattering in pitch angle as shown 
in Figj^t, where an average pitch angle of electrons with respect to Bq and its standard 
deviation are plotted as a function of time. The standard deviation just after its rapid 
growth at Qot ~ 70 is ~ 0.5, which roughly coincides with an analytical estimate of the 
maximum pitch angle width of an electron resonating with the monochromatic parent wave 
(see Appendix). After the parent wave disappears, most of the electrons are detrapped by 
the parent wave and start to wander in the phase space. As time passes, a high energy tail 
appears in the electron distribution function (Fig 12^), which gradually approaches a power- 
law spectrum with index a ~ —2.5 {Qot = 316: gray solid line, 1616: black solid line). At 
the same time the wave amplitude spectrum develops the power-law type spectrum also, 
with the index a ~ —2.0, while the wave spectral peak shifts toward lower wavenumbers. 



A behaviour of the most efficiently accelerated electron is shown in Figl3l FigjSK 
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indicates energy time history of the electron. Spatial trajectory of the electron during 
< not < 600 (800 < not < UOO) is plotted in Figl3t(d). A background color scale 
denotes strength of magnetic fluctuations = y^S|~+~Bf. Horizontal arrows in Fig|3b(d) 
and FigOlL correspond to each other denoting time intervals where strong accelerations oc- 
cur. Figl3^(f) shows a snapshot of the wave profile at the time denoted as the dashed line 
in Figj3l3(d). The red, blue, and black lines represent By, B^, and an envelope, respectively. 
The dotted line indicates the position of the electron. It is seen that strong acceleration 
occurs when the electron is trapped in a trough of the magnetic envelope. Such sharp enve- 
lope troughs are temporarily observed in various regions of the system throughout the run. 
Hence, the electron experiences similar acceleration processes several times. In this example 
one may recognize four such acceleration periods. Only the second and the third periods are 
marked (the first (fourth) one is 160 < Qot < 190 {Qot > 1450)). The acceleration occurs 
always perpendicular to Bq as seen in FigOb which shows a trajectory of the electron in 
the u± — u\\ space, where u± and u\\ denote the four-velocities perpendicular and parallel to 
Bq. The red markers indicate positions of the electron in the phase space at the same time 
as Figl3^ (lower marker) and Fig|3]F (upper marker). We checked a hundred most efficiently 
accelerated electrons' trajectories and confirmed that all of them show essentially the same 
features as mentioned above, i.e., trapping within the envelope troughs and successive per- 
pendicular acceleration. In the next section the acceleration process is modeled and analyzed 
in detail. 



3. Acceleration of High Energy Electrons 

3.1. Coherent waves observed in the simulation: modeling 

In the PIC simulation shown above sharp magnetic envelope troughs are locally formed 
throughout the period of strong electron accelerations. In a successive decay process a num- 
ber of daughter Alfven waves propagating both parallel and antiparallel to Bq are excited. 
Depending on their phases, amplitudes of some wave modes are sometimes locally canceled 
out each other, or they are simply less dominant in amplitude. Then, there appear some 
regions where two oppositely propagating waves dominate. Here, the envelope structures 
observed in the PIC simulation are modeled by a superposition of such two oppositely prop- 
agating waves as follows. 

B,,,coskx( ] 



B^ J \ sinut 
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Eq.([T]) is equivallent to 

g Byj / cos(A;x — ut) \ ^ / cos(— /cx — ut) \ 
^ 2 \ —sm.{kx — cot) J 2 \ — sin(— A;x — cut) J 

Hereafter, we assume u and k are both positive without losing generahty. Typical waveforms 
at two different time domains are shown in FigH^ and b. The black solid and dashed lines 
denote y and z components of rotating carrier waves, and the gray lines represent envelopes 
which is independent of time. In the neighborhood of a trough, it is confirmed that eqs.([l]) 
and ([2]) give a reasonable model of a waveform observed in the PIC simulation shown in 
FigHt. In the following motion of a test particle in this system is analyzed. 



3.2. Motion of a particle in given EM fields 

In this section we consider motion of an electron in the electromagnetic waves given by 
eqs.([T]) and ([2]). The equation of motion is 

(^u e u ^\ , , 

E+-XB , (3) 



dt uiqc \ 7 

dx Ux / . N 

-77 = — C, (4) 

dt 7 

where u = 7v/c, 7 = ^/l + vP, e and mo indicate charge and rest mass of the electron, B = 
-Bqx+Bw, and E = E^, respectively. When we write u = {u\\{t), u_i{t) cos (t>(t),u_i_(t) sin0(t)) 
and introduce normalized variables as ^ = xQq/c, t = Q^t, n = kc/Qo, v = uj/Qq, = 
Bw/Bq, and Vph = i^/n, eqs.(l3]) and (jlj) are written as follows. 

u\\ = by^— cos kC, sin ijj (5) 

7 

w± = &to sin K,^ cos ip cos sin ip^ , (6) 

ip = — — ( sin sin ijj -\ — - cos cos \ — { v ) . (7) 

ul\ 7 J \ 1/ 

e = J. 8 

7 

Here, ip = (p — z/r, and the dot denotes time derivative in terms of r, respectively. Note that 
a variation of normalized particle energy, or the Lorentz factor, is given as 

7 = bujVph— sin cosip. (9) 

7 

In the following behaviours of the electron is discussed by using eqs.([5])-([9]). 



- 7- 



3.3. Fixed point analysis 

The above set of equations clearly has several fixed points as listed in Table [H Here, uq 
denotes a value of u± at a corresponding fixed point and 70 = -y^l + Mq- Apparently, the fixed 
points I - III (IV - VI) correspond to troughs (crests) of the magnetic envelope. As far as the 
acceleration is concerned, it is easily inferred that the fixed points IV - VI are not important 
since the electric field strength is very weak around there. Actually, efficient acceleration 
observed in the PIC simulation always occurs around the troughs of the magnetic envelope. 
Therefore, only the fixed points I - III are focused here. Stability of these fixed points is 
discussed in the following. 



3.3.1. Stability of the fixed points I - III 

Expanding eqs.([5])-(j9]) around the fixed point I and retaining only the first order terms, 
we obtain 

Su\\ = —byj—nS^, (10) 
7o 

6ui_ = -b^VphSij, (11) 

Sip = ^Vph6u± - ^57, (12) 
■"o 7o 

Su\\ , , 

bi = -^ (13) 
7o 

(^7 = -byj—Vph6i/j = —Su±, (14) 
7o 7o 

where S denotes small first order quantities. 
From eqs. (|TO|) and (|T3|) . we have 



Su\\ = —bw-^i^Su\\. (15) 

7o 

The above expression represents a trapping oscillation with trapping frequency Utrap = 
(bwUoK/'jQY^'^, which is rewritten with the original parameters as utrap = {B^kcuo/ BqVLq'^qY/'^ . 

From eqs. (fTT]) . (IT^ . and ffT^ . on the other hand, we obtain 

Sul = ^"'""^^^ fl - b^Vph-^^ 5u±. (16) 
7o V 4 J 

Here, the first term in the parenthesis (or the second term in the right hand side of eq. (fT2l) ) 
arises from the variation of the Lorentz factor, i.e., due to the relativistic effect. In the 
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non-relativistic limit, therefore, eg. (1161) represents a trapping motion in the perpendicular 
momentum space with frequency ujtrap = b^Vph/uQ = {B^/ Bq){(jj /kc){l/uQ). However, in 
the relativistic case with 6^Wph7o/^o < ^""J- (and 5ip) diverges in time. When uq is small 
enough so that 70 ~ 1, the above inequality is hardly satisfied. In such a case, the fixed 
point is stable. But if mq becomes large and the inequality is satisfied, such an electron gains 
transverse energy while it keeps being trapped in the x— direction. We consider this solution 
later more in detail. 

For the fixed point II, equations corresponding to eqs. (fT5l) and (|T6l) are obtained by for- 
mally changing the sign of hw Therefore, the system is unstable for parallel fluctuations while 
it is stable for perpendicular fluctuations. This fixed point is actually conjugate to the rela- 
tivistic fixed point discussed above. Another fixed point (III), conjugate to the nonrelativistic 
one, which should be a saddle point in m_l — t/^ phase space, appears at (m±,'?/') = (0,n7r) as 
singular points of du±/di/j, where n = 0, ±1, ±2, ■ ■ •. 



3.3.2. Features in u\^ — ip phase space 

In order to study perpendicular dynamics let us first consider a reduced system in which 
parallel quantities are fixed at the fixed point, i.e., u\\ = and = n/2. Then we only have 
to consider the following two equations. 

p = 2y/pbuiVph cos ip (17) 
^ = -^Vph sin ip - [u - 



Here, p = u']_ has been introduced. This system has a Hamiltonian defined as 



H{p,i)) = 2^/pb^VphSimlj + z/p - 2a/1 +p, (19) 

where p = OH /dip and ip = —OH /dp are satisfied. Contours of the Hamiltonian for = 1.0, 
Vph = 0.17, and v = 0.11 are represented as gray lines of the upper (linear scale) and lower 
panels (logarithmic scale) in FigJS^. The above parameters are chosen so that the waves 
interacting with the accelerated electron in the second acceleration stage observed in the 
PIC simulation (FigJH^) is appropriately reproduced. It is easily confirmed that a center at 
Ip = tt/2, which is clearly recognized in the lower panel, corresponds to the nonrelativistic 
stable fixed point of the fixed point I and another center at ip = 3it/2, which can be seen 
both in the upper and lower panels, to the relativistic fixed point II discussed above. 

The nonrelativistic center should satisfy the following condition from Table [H 

-^ = .-1 (20) 
uo 7o 
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Although — 1/70 = in the small amplitude limit, it is never satisfied in the nonrelativistic 
case since u < 1. Therefore, the nonrelativistic center appears only when wave amplitude 
becomes finite. The closed trajec tories around this center ess e ntially coincides with the 



nonresonant trapping discussed by lKuramitsu fc Krasnoselskikhl (120051 ) 



At the relativistic center, the resonance condition i/ — I/70 = should be satisfied in 
the ultrarelativistic limit, since b^Vph/uQ is negligible. In such a case relativistic decrease of 
the gyro frequency allows another resonance with low frequency waves. This resonance can 
also be present in small amplitude limit. Relativistic linear resonance conditions between an 
electron and two oppositely propagating waves (7Z/ =f ku\\ — 1 = 0) are shown in FigjSt. The 
two resonance lines intersect at u\\ = where z/ — I/7 = and 7 = yT+"u^ are satisfied, 
while there never appears such an intersection in nonrelativistic limit. This indicates that 
an electron can resonate simultaneously with two waves in the relativistic case. By using 
eq. (fT9|) . the maximum width of the separatrix in u± is estimated as 

Au^ = 4./5. (21) 

V Kl/ 

Furthermore, the maximum u± on the separatrix is given by 





For instance, u±^max ~ 18.4 for the parameters used in FiglHl This gives a good agreement 
with the electron energy achieved in the end of the second acceleration stage observed in the 
PIC simulation. 

So far we have only looked at a subset of the system in which Su\\ = 0. If we allow 6u\\ 
to slightly deviated from zero, it should satisfy 

6u\\ = — by^-^K sin ibSu\\ (23) 
7o 

around the fixed point in the original system eqs.([2])-([S])- Hence, the system is stable for 
< ip < n and is unstable for vr < < 27r in terms of parallel fiuctuations. The black solid 
line in Figl5^ and[5)D shows a numerical solution of eqs.([5])-(j8]). A trajectory of the electron 
initially positioned near the separatrix at ip = (indecated by a small arrow in FigJS^) is 
represented. The electron moves alomost along the separatrix when it is in < ■?/'(< vr). But 
when the electron enters in it < ip < 27r, its trajectory starts deviating from the separatrix. 
And finally it is clearly detrapped when it approaches ip = 37r/2. The trajectory in the 
momentum space shown in Fig. [Sjo is very similar to what is observed in the PIC simulation 
which is again plotted as a gray line. 
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3.4. Statistics of high energy particles 



Let us briefly discuss the power-law like distribution of electorns in the current accel- 
eration process. It is shown in Figl2^ that a high energy tail with power-law index ~ —2.5 
evolves in the late stage of the simulation. This power-law nature may be related with 
power-law spectrum of wave amplitude of magnetic fluctuations. In the stage of successive 
decay instabilities {flot > 100) a spectral cascade of the magnetic fluctuations occurs with 
keeping a power-law index of the wave amplitudes almost constant with ~ —2.0 (Fig|2]D). 
The enhancement of the wave amplitude at higher wavenumber regime, kc/Qo > 3, is the 
remnant of initial parent wave and side band waves generated by the modulational instabil- 
ities in the early stage of the run. In order to confirm correlations between the wave and 
electron energy spectra, the following test particle simulation with periodic boundary con- 
ditions is performed. Waves are given by superposition of a number of right hand polarized 
Alfven waves with power-law spectrum shown in Figj6] (power-law index a is an external 
parameter), which models the lower wavenumber part of Figl2)D. Both positive and negative 
wavenumber modes are evenly distributed. The wave form is given by 



where Vph/c = 0.17 and {$fc} are initial wave phases which are randomly distributed at 
t = 0, and corresponding transverse electric fields are given by k x E = {\k\vph/c)'B. The 
power-law index a is fixed to 1.0 for kmin k < Kq, where kmin = 27r/L and Kq is the 
coherence wavenumber. The system size L = SOttc/Qq is common for all the following runs, 
and = 0.75 and Kqc/Qq = 0.1 if not specified. 

In such a system where a number of waves propagate in two opposite directions, particles 
are expected to be stochastically accelerated via second order Fermi like process. First, 
this is confirmed by solving nonrelativistic equations of motion of 10^ electrons. An initial 
distribution function is given as a spatially homogeneous gyrotropic ring distribution with 
vj_/vph = 0.1. FiglT] shows distributions in (a) v± — v\\ and (b) X — Wkin phase spaces, and 
(c) energy distribution functions at two different times, Q^t = 1000 and 3000, respectively, 
where v = (fy, v±) denotes particle velocity and Wkin = 'y^/2fp/i is normalized kinetic energy. 
The solid gray lines in FigjTb denote envelope profiles of given magnetic fluctuations. The 
particles are spatially bunched around the troughs of the magnetic envelopes because of 
mirror effect. In the velocity space the particles rapidly pitch angle diffuse (FigJT^.) while 
the averaged particle energy slowly increases in time (FigjTt). These are the features of 
a stochastic or a second order Fermi acceleration. In comparison, drastic changes occur 
if relativistic effects are taken into account. The simulation results with the same initial 




(24) 



- 11 - 



conditions as in the nonrelativistic case are plotted in Figj8]with the same format as FigJTl 
Here, the normahzed kinetic energy is defined as Wkm = {l — Compared with the 

nonlerativistic case, the maximum particle energies at each corresponding time are extremely 
higher, and they are at the same order as that obtained from eq. (l22|) with k = Kq and 
bu, = 0.75 (the black solid and black dashed lines in FigjHt). Spatial bunching of particles 
similar to the nonrelativistic case is found in FigjHjD, while a fraction of the bunched particles 
are accelerated to extremely high energies. Such high energy particles have rather large pitch 
angles, denoting perpendicular acceleration (FiglH^). Interestingly, a power-law like energy 
distribution appears only at QqI = 1000 in FiglHt, although the corresponding energy range 
in the distribution function is small (the black solid line). The power-law index does not 
change even in the case with a different initial ring velocity, v±/vph = 5.1 (not shown). 
However, the spectrum becomes softer when a = —3.0 is assumed (the gray dashed line 
in FiglHt). When the coherence wavenumber is large as Kqc/Qq = 0.4, the bulk electrons 
are accelerated and the energy distribution is no longer power-law at flot = 1000 (the gray 
dotted line in FigJHt). It is noted from eqs. (!2Tl) and (122|) that the minimum u± of the 
separatrix of the relativistic resonance of the dominant wave mode at A; = Kq decreases with 
increasing Kq. Hence, most of electrons which initially distribute below the separatrix can 
enter inside the separatrix through stochastic motions at rather early stage, and they can 
be perpendicularly accelerated within a short time in the similar way seen in the previous 
section (cf. Figj5]). The peak energy is roughly consistent with the u±,maa:(~ 940) estimated 
by using = 0.75 and k = 0.4 from eq. fl22|) . In contrast, for small wave amplitude, b^ = 0.1, 
the maximum energy at Qot = 1000 is much smaller than the value obtained from eq.( l22i) 
(not shown) . The reason may be that in this run the minimum u± of the separatrix is rather 
high because of its narrow width so that particles which initially distribute far below the 
separatrix have not entered in it until this time. 

In FiglHt the power-law like nature appears as a transient state in the system at 
flot = 1000 for a = —2.0. In this run, at later time {flot = 3000), high energy end of 
the particle distribution is so enhanced that the energy distribution does not fit the power- 
law spectrum (black dashed line). It is confirmed that the hump of the high energy part 
grows at least till Q^t = 6000. After sufficiently long time, it probably results in bulk accel- 
eration as seen in the case of small Kq. In the PIC simulation the wave spectrum is not time 
stationary but cascading through the successive decay instabilities as already mentioned. In 
other words, a wave with a certain wavenumber has finite life time. This may be why the 
high energy tail evolves without extra accumulation at high energy end of the distribution 
function in the late stage, Qot > 316, of the PIC simulation (FigJ2t). To confirm this, an ad- 
ditional test particle simulation is performed by assuming that wavenumber of the magnetic 
field with maximum intensity varies in time as Kqc/Qq ~ (f2o^) which gives a reasonable 
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fit in the late stage of tlie PIC simulation. The particle distribution and waves at ^IqI = 1000 
in the run corresponding to the black solid line in FigJHt are chosen as initial conditions. 
Then, in the later time flat = 3000 high energy tail extends roughly obeying the power-law 
and some low energy particles remain unaccelerated, as shown in FiglHt as the gray solid 
line. 



4. Summary and Discussions 

In the present paper an efficient particle acceleration process in the course of successive 
parametric instabilities of large amplitude Alfven waves was investigated. The acceleration 
takes place as a result of interactions between coherent waves in the developing Alfven 
turbulence and relativistic particles. An important point is that relativistic wave-particle 
interactions allow simultaneous resonance between a particle and two different waves. The 
maximum attainable energy through this acceleration process was analytically estimated. 

In this acceleration process a high energy particle is preferentially accelerated. Dur- 
ing the successive decay instabilities, a peak of the wave Fourier spectrum shifts in time 
toward a lower frequency (longer wavelength) regime. Because of relativistic decrease of 
particle's gyro frequency, low frequency waves preferentially resonate with and accelerate 
particles with large energy. Therefore, if once a particle is accelerated and its effective mass 
is increased, in later time such a particle easilt satisfies the resonance condition with lower 
frequency waves generated by successive decay processes. Fig|3^ shows an example of such 
an electron's energy time history in which one can recognize four acceleration phases around 
riot ~ 170,430, 1200, and 1550. 

By utilizing test particle simulations in section 4, it is shown that relativistic effetcs are 
essential in producing nonthermal pa rticles. Note t hat th e system discussed in section 3 is 
exactly the same as that discussed by IWykes et al. J2001al jbl). although they mainly focused 



on stochastic off-resonance diffusion in the phase space with an applicat ion to planetary mag - 



netospheres in mind. Actually, the regular trajectories in Fig.3(d) in IWykes et al.l (l2001al ) 
correspond to the relativistic resonance focused in this paper. Enhancement of Lyapunov 
exponent at ~ Qo/wy, where n is an integer, shown in their Fig.8 might be related with the 
relativistic resonance. In our test particle simulation the power-law like energy distribution 
function similar to what was observed in the PIC simulation is reproduced by assuming the 
time evolving Fourier spectrum of Alfven turbulence. Finite life time of a wave mode in 
developing turbulence may contribute to creation of such a distribution function. However, 
details of the acceleration process including analytical estimate of the power-law index are 
still unclear and will be investigated in the near future. 
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The acceleration process discussed in this paper is essentially different from some of the 
recently studied coherent acceleration processes, i.e., electron surfing acceleration and wake- 
field acceleration. In these processes electrostatic field plays essential roles. Furthermore, 
the processes mainly affect electrons since generated electrostatic waves have rather high 
frequencies, while ion accel eration may also occur in a very late stage of nonlinear evolution 
of a system (lHoshindl2008l ). In the acceleration process discussed here roles of electrostatic 
fields are subdominant, although they are necessary for the decay instabilities. It should be 
noted further that the process may be able to act on ions too when a left hand polarized 
wave is introduced as an initial parent wave. The Alfven waves are essentially incompressible 
so that they may survive for rather long period compared with electrostatic Langmuir or ion 
acoustic waves. Hence, the process may last for long time, and may also follow the above 
mentioned electrostatic acceleration processes in some occasions. 

In the present study all the multidimensional effects have been excluded. In higher di- 
mensional cases daughter waves propagating in ob lique to the ambient magnetic field can also 
have finite growth rates (jVinas fc Goldsteinlll99ll ). It is confirmed by MHD sir nulation that 



these waves destroy the planar stru cture assured in a one dimensioal simulation (iGhosh et al. 



19931 . 1 1994 iDel Zanna et al.ll200ll ). However, in low frequenc y regime the decay instability 
of parallel progagation is dominant (jVinas &: Goldsteinlll99ll ). Therefore, at least in such a 
regime the acceleration process observed here are expected to work, while acceleration rate 
may decrease to some degree because of appearance of nonplanar structures. This is similar 
to what is discussed for electron surfing acceleration in which th e acceleration takes plac e 
even in two dimensional cases despite decrease of acceleration rate (lAmano fc Hoshindl2008l ) . 
At all events, further investigations are necessary for estimate efficiency of this acceleration 
process in a more realistic situation. 

We finally give a comment on applications of this acceleration process. Since the situ- 
ation simulated in section 2 is rather general, there may be several fields of application like 
a pulsar wind nebula, a pulsar magnetosphere, an outflow of GRB, and an AGN jet, and 
so on. Vicinity of a relativistic shock is probably one of the candidates. Similary to the 
earth's foreshock, large amplitude Alfven waves may be generated by beam-plasma interac- 
tions between backstreaming ions and an upstream plasma, and the waves nonlineary evolve 
via parametric instabilities. Since the acceleration is efficient and locally takes place, it may 
contribute to the so-called injection process into the DSA. As another possible case, cosmic 
ray-plasma interactions upstrear n of a supernova remnant shock are now exten sively stud- 
ied a fter the pioneering works bv llucek fc Bel3 J200oh . lBell fc Lucekl J200lh . andEsS fl2004j . 
20051 ). Most of such studies pay attention to amplification of upstream magnetic fluctuations 
that can be scatterers of cosmic rays. And that is expected to result in increase of maximum 
attainable energy in the DSA process. On the other hand, amplified magnetic fluctuations 
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may lead to the local and the coherent acceleration process discussed here. In other words, 
some particles may be accelerated through this process without crossing a shock. In this 
sense the process is similar to the so-called second order Fermi acceleration, although an 
efficiency of the acceleration process discussed here is much higher than that of the second 
order Fermi process as shown by the test particle simulation. 

The authors thank Victor Muiioz for useful discussions. The PIC simulation was per- 
formed by the super computer in ISAS/JAXA Sagamihara. This work was supported in part 
by Incentive aid to the prominent research. Interdisciplinary graduate school of engineering 
sciences, Kyushu University 2007. 



A. Motion of a relativistic particle in a monochromatic wave 



According to iKuramitsu &: Krasnoselskikhl (120051 ) , an equation of motion of an electron 
in a monochromatic circularly polarized wave in a wave frame is reduced as 

dH ^_dH 



(Al) 



where /i is the pitch angle cosine, ip the particle gyrophase with respect to the wave phase, 
and 

K / 1 \ ^ 

i/ = yU+-j +6„Vl-/i2 cosV- (A2) 

is Hamiltonian. The above expressions are held even when relativistic effects are taken 
into account by putting K = ukc/Qo and by^ = {B^/ Bq) ^1 — f^^/ c^. Electron trajectries 
interacting with the parent wave {B^/Bq = 1 and kc/Qo = 6.21) for different values of H 
are plotted in FiglH Here, u = 0.64 is assumed since the value is derived as average one at 
flot = 70 in the PIC simulation. The factor of f^^/c^ is neglected b e cause of its smallness. 
Figiniis essentially the same as Fig. 5 in IKuramitsu fc Krasnoselskikhl (120051 ). The maximum 
value of half the width of the separatrix is ~ 0.5 (at ip = it), which roughly coincides with 
the standard deviation of pitch angle observed in the PIC simulation at QqI = 70 (Figl^t). 
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Table 1: Fixed points 
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Fig. 1. — Time evolution of wave power spectra. The right and left hand helicity modes 
of Bz component are Fourier decomposed in left and middle panels, respectively. The right 
panel shows Ex component. See the electronic edition of the Journal for a color version of 
this figure. 
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Fig. 3. — Behaviours of an accelerated electron, (a) Time history of the electron energy. 

(b) Trajectories of the electron in u± — u\\ space and (c)(d) in X — t space. The arrows in 

(c) and (d) indicate time domains where strong acceleration occur as shown by the arrows 
in (a). Background color scale in (c) and (d) denotes amplitude of magnetic fluctuations. 
(e)(f) Spatial profiles of magnetic fiuctuations in two acceleration phases represented by 
dashed lines in (c) and (d). The red, blue, and black solid lines show By, components 
and envelope, respectively. The dotted fines sfiow positions of tfie particle at correspponding 
times. See tfie electronic edition of tfie Journal for a color version of tfiis figure. 
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Fig. 4. — (a), (b) Modeled wave forms at different time domains (i?,^, = l,u/kc— 0.2), and 
(c) an example of a waveform seen in the PIC simulation at VL^t = 1170. 
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Fig. 5. — Electron trajectories in (a) u^ — il), (b) — phase space, and (c) linear resonance 
conditions of wave-particle interactions in relativistic (solid line) and nonrelativistic (dashed 
line) plamsas. Solutions of eg. (1191) are plotted as gray lines in (a). In (b) a gray line 
denotes a trajectory of the electron shown in Fig|3l A numerical solution of eqs.([5]) - ([8]) for 
n^/f2o = 1, fcc/|f2o| = 0.65, and ct;/|i7o| = 0.11 is indicated as black lines in (a) and (b). 
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Fig. 6. — Wave power spectrum used in the test particle simulation. 
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Fig. 7. — Results of nonrelativistic test particle simulation. Electron distributions in (a) 
v± — v\\ and (b) X — W^n phase spaces, and (c) energy distribution functions at — 1000 
and 3000, respectively. The solid gray lines in (b) denote profiles of magnetic field envelopes 
at the corresponding times. 
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Fig. 8. — Results of relativistic test particle simulation in the same format as Fig 13 A 
number of energy distribution functions are plotted in (c). See datails in the text. 
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Fig. 9. — Trajectories of relativistic electrons interacting with the monochromatic circularly 
polarized parent wave. B^)/ Bq = 1, kc/Qo = 6.21, and u = 0.64 are assumed. 



